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It was recently shown that conventional phonon-electron interactions induce triplet pairing states 
in time-reversal invariant 3D Dirac semi - metals provided magnetic impurities or exchange inter¬ 
actions are strong enough 1 . The order parameter in this case is a vector field. Starting from the 
microscopic model of the isotropic Dirac semi-metal, the Ginzburg-Landau energy for this field is 
derived using the Gor’kov technique. It was found that the transversal coherence length is much 
smaller than the longitudinal, £ L = 4\/2£ T , despite the isotropy. Several new features appear when 
an external field is applied. The Ginzburg - Landau model is used to investigate magnetic properties 
of such superconductors. Using the small deviation method the magnetic penetration depth was 
found also to be vastly different for longitudinal and transverse fluctuations Xt/Xl = 4\[2. As a 
result the superconductor responds as type I to a transverse perturbation, while with respect to 
a longitudinal perturbation it behaves as type II. At large fields the order parameter orients itself 
perpendicular to the field direction. The triplet superconductor persists at arbitrarily high magnetic 
field (no upper critical magnetic field) like in some p wave superconductors. 

PACS numbers: 74.20.Fg, 74.90.+n, 74.20.Op 


I. INTRODUCTION 

Recently 3D Dirac semi - metals (DSM) like Na^Bi 
and CdzAs 2 with electronic states described by Bloch 
wave functions, obeying the ’’pseudo-relativistic” Dirac 
equation (with the Fermi velocity vf replacing the veloc¬ 
ity of light) were observecP and attracted widespread at¬ 
tention. The discovery of the 3D Dirac materials makes 
it possible to study their physics including remarkable 
electronic properties. This is rich in new phenomena 
like giant diamagnetism that diverges logarithmically 
when the chemical potential approaches the 3D Dirac 
point, a linear-in-frequency AC conductivity that has 
an imaginary parfP, quantum magnetoresistance show¬ 
ing linear field dependence in the bulk^l. Most of the 
properties of these new materials were measured at rela¬ 
tively high temperatures. However recent experiments 
at low temperature on topological insulators and sus¬ 
pected 3D Dirac semi-metals exhibit superconductivity. 
Early attempts to either induce or discover supercon¬ 
ductivity in Dirac matrials were promising. The well 
known topological insulator Bi 2 Se^, doped with Cu , be¬ 
comes superconducting at T c = 3.8A®. At present 
its pairing symmetry is unknown. Some experimental 
evidence 6 point to a conventional phononic pairing mech¬ 
anism. The spin independent part of the effective elec¬ 
tron - electron interaction due to phonons was studied 
theoretically '. For a conventional parabolic dispersion re¬ 
lation, typically independent of spin, the phonon mech¬ 
anism leads to the s-wave superconductivity. The lay¬ 
ered, non-centrosymmetric heavy element PbTaSe 2 was 
found to be superconducting. Its electronic proper¬ 
ties like specific heat, electrical resistivity, and magnetic- 
susceptibility indicate that PbTaSe 2 is a moderately cou¬ 
pled, type-II BCS superconductor with large electron- 
phonon coupling constant of A = 0.74. It was shown the¬ 
oretically to possess a very asymmetric 3D Dirac point 
created by strong spin-orbit coupling. If the 3D is con¬ 


firmed, it might indicate that the superconductivity is 
conventional phonon mediated. 

More recently when the Cu doped Bi 2 Se 3 was sub¬ 
jected to pressur^l, T c increased to 7 K at 30 GPa. Quasi- 
linear temperature dependence of the upper critical field 
H c 2 , exceeding the orbital and Pauli limits for the sin¬ 
glet pairing, points to the triplet superconductivity. The 
band structure of the superconducting compounds is ap¬ 
parently not very different from its parent compound 
Bi 2 Se 3 , so that one can keep the two band k • p de¬ 
scription (Se p z orbitals on the top and bottom layer of 
the unit cell mixed with its neighboring Bi p z orbital). 
Electronic-structure calculations and experiments on the 
compounds under pressurd^ reveal a single bulk three- 
dimensional Dirac cone like in Bi with large spin-orbit 
coupling. Moreover very recently some pnictides were 
identified as exhibiting Dirac spectrum®. This effort re¬ 
cently culminated in discovery of superconductivity in 
CdsAs 2 P. It is claimed that the superconductivity is p- 
wave at least on the surface. 

The case of the Dirac semi-metals is very special due 
to the strong spin dependence of the itinerant electrons’ 
effective Hamiltonian. It was pointed out23 that in this 
case the triplet possibility can arise although the triplet 
gap is smaller than that of the singlet, the difference 
sometimes is not large for spin independent electron - 
electron interactions. Very recently the spin dependent 
part of the phonon induced electron - electron interac¬ 
tion was considered 11 and it was shown that the singlet 
is still favored over the triplet pairing. Another essential 
spin dependent effective electron-electron interaction is 
the Stoner exchange among itinerant electrons leading to 
ferromagnetism in transition metals. While in the best 
3D Weyl semi-metal candidates it is too small to form a 
ferromagnetic state, it might be important to determine 
the nature of the superconducting condensate. It turns 
out that it favors the triplet pairing 12 . Also a modest 
concentration of magnetic impurities makes the triplet 
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ground state stable. 

As mentioned above generally the applied magnetic 
field is an ultimate technique to probe the supercon¬ 
ducting state. In a growing number of experiments, 
in addition to magnetotransport, magnetization curves, 
the magnetic penetration depth and upper critical mag¬ 
netic field were measured 1 ®. It is therefore of importance 
to construct a Ginzburg - Landau (GL) description — 1 
of these novel materials to study inhomogeneous order 
parameter configurations (junctions, boundaries, etc.) 
and magnetic properties that typically involve inhomo¬ 
geneous configurations (like vortices) not amenable to a 
microscopic description. 

In the present paper we derive such a GL type theory 
for triplet superconductor from the microscopic isotropic 
DSM model with attractive local interaction. The order 
parameter in this case is a vector field and the GL theory 
of vector field already considered in literatur^^ 17i in con¬ 
nection with putative p - wave superconductors have sev¬ 
eral extraordinary features, both quantitative and quali¬ 
tative. 

The paper is organized as follows. The model of the 
(phonon mediated or unconventional) local interactions 
of 3D Dirac fermion is presented and the method of its 
solution (in the Gorkov equations form) including the 
symmetry analysis of possible pairing channels and the 
vectorial nature of the triplet order parameter is given 
in Section II. In Section III the Gorkov formalism, suf¬ 
ficiently general to derive the GL equations, is briefly 
presented. The most general form of the GL energy of 
the triplet superconductor in magnetic held consistent 
with the symmetries is given in IV. The coefficient of the 
relevant terms are calculated from the microscopic DSM 
model in section V. Section VI is devoted to applications 
of the GL model. The ground state degeneracy, the char¬ 
acter of its excitations and basic magnetic properties are 
discussed. The vector order parameter is akin to optical 
phonons with sharp distinction between transverse and 
longitudinal modess. Transverse and longitudinal coher¬ 
ence lengths and penetration depths are calculated and 
the upper critical magnetic held is discussed. Section VI 
includes generalizations to include Pauli paramagnetism, 
discussion of an experimental possibility of observation 
of the excitation and conclusion. 


II. THE LOCAL PAIRING MODEL IN THE 
DIRAC SEMI-METAL. 

A. Pairing Hamiltonian in the Dirac semi-metal. 

Electrons in the 3D Dirac semimetal are described 
by held operators ipf s (r), where f = L,R are the val¬ 
ley index (pseudospin) for the left/right chirality bands 
with spin projections taking the values s =f, i with re¬ 
spect to, for example, the 2 axis. To use the Dirac 
(’’pseudo-relativistic”) notations, these are combined into 
a four component bi-spinor creation operator, ip' = 


V 4 j.) > wh ose index 7 = {/, s} takes four 
values. The non-interacting massless Hamiltonian with 
Fermi velocity vf and chemical potential p readd 1 ® 

K = J ip + (r) Kip (r); (1) 

K^yS — tfVUp' V QLry§ flSsyS , (2) 

where the three 4x4 matrices, i = x, ?/, z, 

<*> 

are presented in the block form via Pauli matrices u. 
They are related to the Dirac 7 matrices (in the chiral 
representation, sometimes termed ’’spinor”) by a = 
with 

*-(SS)- (4) 

Here 1 is the 2x2 identity matrix. 

We consider a special case of 3D rotational symmetry 
that in particular has an isotropic Fermi velocity. More¬ 
over we assume time reversal, Qip (r) = ia y ip* (r), and 
inversion symmetries although the pseudo- Lorentz sym¬ 
metry will be explicitly broken by interactions. The spec¬ 
trum of single particle excitations is linear. The chemical 
potential p is counted from the Dirac point. 

As usual in certain cases the actual interaction can be 
approximated by a model local one: 

Veff = - 1 J^air) (r) ipp (r) ip a (r). (5) 

Unlike the free Hamiltonian K , Eq. ([!]), this interaction 
Hamiltonian does not mix different spin components. 
Spin density in Dirac semi-metal has the form 

s M = \tp + (r) Vip (r), (6) 

where the matrices 

s = - Q rs = ( 0 Z ) ' < 7 > 

75 = ( 0 1) - 

are also the rotation generators. 

B. The symmetry classification of possible pairing 
channels. 

Since we consider the local interactions as dominant, 
the superconducting condensate (the off-diagonal order 
parameter) will be local 

O = J ipt (0 Mapip'p (r), (9) 
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where the constant matrix M should be a 4 x 4 anti¬ 
symmetric matrix. Due to the rotation symmetry they 
transform covariantly under infinitesimal rotations gen¬ 
erated by the spin S % operator, Eq.([ 6 |: 


/ ( r ) M ocf}i>p (r), 0 7 ( r ') (r') ( 10 ) 

J r,r' J 

= - / 0 + (r) m s M SK + ) 0 + (r). ( 11 ) 

J r 


obey the Gor’kov equationd^ 


_ dG,«^r-y.T 0 _£ _, r |^|^ Gjj>(r „, r;rV) 

(17) 

( r > u r > r ) F £ K 0>u r ', r ') = <^ 7K ^ ( r - r ') ^ ( r - r '); 

9F+ k (r,r;r ' ,t') 

8t 


J ( r |^ 7/31 r ") F /£c ( r ",r;rV) 


(r, r; r, r) G^ (r, r, r', r') = 0 . 


Here and in what follows ”t” denotes the transpose ma¬ 
trix. The representations of the rotation group therefore 
characterize various possible superconducting phases. 

Out of 16 matrices of the four dimensional Clifford 
algebra six are antisymmetric and one finds one vector 
and three scalar multiplets of the rotation group. The 
multiplets contain: 

(i) a triplet of order parameters: 


These equations are conveniently presented in matrix 
form (superscript t denotes transposed and I - the iden¬ 
tity matrix): 



(X, X") G ( X ", X') - A (A') F+ (X, X') = IS (X - 



(X, X") F+ (X", X') + A 4 * (A) G {X, X') = 0. 


= {-/3a z ,-if3'y 5 ,/3a x } = {T x ,T y ,T z } 
The algebra is 

+ TjYj\ = 2 iSijkTk- 


Note that the three matrices Ti are Hermitian. 
(ii) three singlets 


( 12 ) 

(13) 


(14) 


Here X = (r, r), A a0 (X) = gF 0a (X, X) and 


D- 1 0 (X,X') = -S a0 -S(X 


X')—S (t 


T ') 




(19) 

In the homogeneous case the Gor’kov equations for 
Fourier components of the Green’s functions simplify 
considerably: 


Mf = ia y ; M (f = i E y ; = —i/3ayj 5 . (15) 

Which one of the condensates is realized at zero temper¬ 
ature is determined by the parameters of the Hamilto¬ 
nian and is addressed next within the Gaussian approx¬ 
imation. As was shown in our previous worll 1 ^ either 
exchange interactions or magnetic impurities make the 
triplet state a leading superconducting channel in these 
materials. Therefore we will consider in the next section 
only the vector channel. 


III. GORKOV EQUATIONS AND THE 
TRIPLET PAIRING 

A. Gorkov equations for Green’s functions in 
matrix form 

Using the standard BCS formalism, the Matsubara 
Green’s functions (r is the Matsubara time) 

Gat 3 ( r > r ; r', r) = - (r T ip a (r, r) ^ (r', r')) ; (16) 
Fat 3 (U u r '> T ') = ( T ri>a ( r , T ) 0/3 ( r ': T ')) 5 
( r > F r', r') = (r T ipl (r, r) (r 4 , r')) , 


D 1 (w, p) G (w, p) — AF + (w, p) = /; (20) 

D- 1 (u, p) F + (u, p) + A 4 *G (u>, p) = 0. 

The matrix gap function can be chosen as 

A/ 3 7 = gFj 0 (0) = A Z M J0 , (21) 

with real constant A z . Here D _1 (w, p) = ioj + g — a ■ p, 
is the noninteracting inverse Dirac Green’s function for 
the Hamiltonian Eq.|l]) and D~ x (u, p) = iui — g — a : 4 p, 
where u> n = irT (2n + 1) is the fermionic Matsubara fre¬ 
quency. 

Solving these equations one obtains (in matrix form) 

G -1 = TV 1 + ADA 4 *; (22) 

F+ = -DA 4 *G, 

with the gap function to be found from the consistency 
condition 

A 4 * = - j^i}A 4 *G. (23) 

cop 

Now we find solutions of this equation for each of the 
possible superconducting phases. 
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B. Homogeneous triplet solution of the gap 
equation. 


A. The GL description for a vector order 
parameter 


In this phase rotational symmetry is spontaneously 
broken simultaneously with the electric charge U (1) 
(global gauge invariance) symmetry. Assuming z direc¬ 
tion of the p - wave condensate the order parameter ma¬ 
trix takes a form: 


A = A Z T Z = A z /3 a x 


(24) 


where A z is a constant. The energy scale will be set by 
the Debye cutoff To of the electron - phonon interactions, 
see below. 

The spectrum of elementary excitations at zero tem¬ 
perature was discussed in refP^. There is a saddle point 
with energy gap 2A Z on the circle p\ = p 2 +p 2 = g 2 /v F , 
p z = 0. The gap A z as a function of the dimensionless 
phonon-electron coupling A = gN , where N being the 
density of states (all spins and valleys included), increases 
upon reduction in g. At large g » To, as in BCS, the 
gap becomes independent of g and one has the relation 

(25 > 

leading to an exponential gap dependence on A when it 
is small: A z = To /sinh (12/A) ~ 2Toe~ 12 ^ x . 


The critical temperature is obtained from Eq. (23) with 


disctretized w by substituting A z = 0. To utilize the 
orthonormality of Ti, Tr (TiZ)*) = 4 Sij, one multiplies 
the gap equation by the matrix T z /g and takes the trace: 


-Tr {T Z T*) = 4 = T C B Z 
9 9 


The bubble integral is 


Bij Y. Tr ( TiDT*D) - 4 SijTr x 


cjp 


E 


V F ( P± ~ Pi ) + M 2 + 


<*4 + (v 2 F p 2 - g 2 ) +2c4 (v 2 F p 2 + g 2 ) 


(26) 


(27) 


{.28) 


Performing first the sum over Matsubara frequencies and 
then integrate over q one obtains, similarly to the singlet 
BCS, (see Appendix A for details): 


^lEr 
7T 


T c = alslT D e- 12 ' x , (29) 

where logy^ = 0.577 is the Euler constant. 


The static phenomenological description is determined 
by the GL free energy functional F [A (r) , A (r)] ex¬ 
panded to second order in gradients and fourth order 
in A. In a magnetic field B, as usual, space derivatives 
of the microscopic Hamiltonian become covariant deriva¬ 
tives V V+ifeA, e* = 2e due to gauge invariance 

under Aj —>■ e lx ^ T ' 1 A\,Ai —> Ai — ffVx- Naively the only 
modification of the GL energy is in the gradient term, 
Eq. (30); the most general gradient term consistent with 


rotation symmetry and the U (1) gauge symmetry is 


Agrad = A 


u T {(DjAi)* (DjAi) - (DiAj)* (DjAi)} 
+u L (DiAj)* (DjAi) 

(30) 

The factor N, the density of states, is customarily in¬ 
troduced into energjJni. It was noted iiP, that, unlike 
in the usual scalar order parameter case, the longitudi¬ 
nal and transverse coefficients are in general different, 
leading to two distinct coherence lengths, see Section V. 
Possibilities for the local terms are^ 


Fioc = N (M) ^ {a (T - T c ) A*Ai + ^ (A* A if + ^ |A»Ai 

( 31 ) 

The magnetic part, F mag = B 2 /8n, completes the free 
energy. 


B. The set of GL equations for triplet order 
parameter 

The set of the GL equations corresponding to this en¬ 
ergy are obtained by variation with respect to 4/* and 
Ai. The first is: 


l u t ( $ijD 0 {Di,Djj ) + 0 ul {Di,Dj} ^ Aj + 

(32) 


+a (T - T c ) Ai + PiAiAjAj + P 2 A*AjAj = 0. 

The anticom mu tator appears due complex conjugate 
terms in Eq.(30j 20 . The Maxwell equation for the su¬ 


percurrent density is: 


IV. A GENERAL GL DESCRIPTION OF A 
TRIPLET SUPERCONDUCTOR IN A 
MAGNETIC FIELD. 


Ji = -j—N (g) {u T A*DiAj + iiA]D ; A. } + cc, (33) 


In this section the effective description of the super¬ 
conducting condensate in terms of the varying (on the 
mesoscopic scale) order complex parameter vector field 
A i (r) is presented. 


where u = Ul — ut- 

Having established coefficients ut,l,Pi 2 > Pz an d a, 
our aim in the next Section is to deduce them from the 
microscopic Dirac senri-nretal model. 
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V. GL COEFFICIENTS FROM THE GOR’KOV 
EQUATIONS 


For the calculation of coefficients of the local part, the 
homogeneous Gor’kov equation, Eq. (23) suffices, while 
for calculation of the gradient terms a general linearized 
equation, Eq. (18) is necessary. Magnetic field effects are 
required only for the calculation of the Zeeman term co¬ 
efficients and normalization of the order parameter via 
supercurrent density. 


A. Local (potential) terms in Gor’kov 


Iterating once the equation Eq. (23) with help of 
Eq. (22) one obtains the local terms to third order in the 
gap function: 


In particular, 


-A * 4 (X) = D+(X,X) (43) 

9 



D 4 (X - X') A** (X 1 ) D (X 1 


A()14) 


The anomalous Green’s functions are no longer space 
translation invariant, so that the following Fourier trans¬ 
form is required: The (time independent) order parame¬ 
ter is also represented via Fourier components A* (X) = 
e _lP ' r A* (P). The linear part Gor’kov equation (this 
time including nonlocal parts) clS the only ’’external” 
momentum P reads 


-A* 4 (P) + E D (w, p) A 4 * (P) D (w, p - P) . (45) 

9 


-A*‘+V I ~ £) (w,p)A* 4 Djw,p)- 

9 ^ \ -D (w, p) A* 4 D (w, p) AD (w, p) A*‘D (w, p) 

(34) 

Using A 4 * = A*Tj, multiplying by T 4 and taking the 
trace, one gets the linear local terms 


Na (T — T c ) A* = -A*, 
9 


(35) 


where the bubble integral was given in Eq.(27). Express¬ 


ing g via T c , see Eq. (|29j) allows to write the coefficient a 
in the Gork 

a(T — T c ) = 


of A* in the Gorkov equation Eq. (23) is 
8^ 2 


T 


3Tr 2 v%h 3 N g T, 


4T~T C 
3 T c 


(36) 


where N is the density of states. 

The cubic terms in Eq.(34), multiplied again by T 4 and 
’’traced” take the form 

(37) 


N (PiAjAjA* +/3 2 A*A*A,;) 


= -A*A fc A? 


2>{ 


TjDTf* DT k DTf* D 


}. (38) 


The calculation is given in Appendix A and results in: 


01 = 


7C(3) 1 


^ = -3^ 


20t r 2 T 2 

The Riemann zeta function is £ (3) = 1. 2 . 


(39) 


To find the coefficients of the gradient terms, one 
should consider contributions quadratic in P from the 
expansion of both A 4 * (P) a nd D (oj, p — P). In view of 
the gap equation Eq.( 26|27 l, the expansions of A 4 * (P) 
cancel each other up to small corrections of order T — T c . 
So that multiplying by T\ and taking the trace 


\p k Pi E Tr { T ' D p) T T D 'ki ( w - p)} 


(46) 


where D kl (w,p) = 9 ■ Comparing this with the 
gradient terms in the GL equation, Eq.(32), see Appendix 
B for details, one deduces 


Ut 


28C (3 )v 2 F h 2 _ 1 
15t r 2 Tf’ 


(47) 


Note the very small longitudinal coefficient, << ut- 
As we shall see in the following section it has profound 
phenomenological consequences. 


VI. BASIC PROPERTIES OF THE TRIPLET 
SUPERCONDUCTOR 


A. Ground state structure and degeneracy 


B. Linear gradient terms 


To calculate the gradient terms, one first linearizes the 
Gor’kov equations, Eq.(18) 


[ D -1 (X, X") G (X", X') 
Jx" 

= 16 (X - X') -+G = D~ 1 ; 

F+ (X, X 4 ) 


(40) 

(41) 


[ D 4 (X - X") A* 4 (X") D (X" - X'). (42) 
J X" 


A ground state is characterized by three independent 
parameters corresponding to three Goldstone bosons. 
The GL energy is invariant under both the vector SO( 3) 
rotations, A, —> RijXj , and the superconducting phase 
1/(1), Aj -A e* x A;. In the superconducting state char¬ 
acterized by the vector order parameter A (|A| = A, 
energy gap) the U (1) is broken: U (1) -A 1, while the 
SO(3) is only partially broken down to its SO(2). There 
are therefore three Goldstone modes. Here we explicitly 
parametrize these degrees of freedom by phases following 
refP. Generally a complex vector field can be written as 
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A = A (n cos (j)+im sin </>), 


(48) 


where n and m are arbitrary unit vectors and 0 < <j> < 
7t/2. Using this parametrization the homogeneous part 
of the free-energy density, Eq.(31), takes the form 


For a perturbation with wave vector k the energy to 
quadratic order in fluctuations reads: 


SF - ST' 

N ~ 


a £+£ } {2« (T c - T) e* £ + u T W T + uW} 


where 


8ne* 2 N 


A* ( k 2 Sij - kikj) Aj 


(54) 


floe = N 


a(T — T c ) A 2 

+ \ (/?i + 0 2 ( cos2 (20) + ( n ' m ) 2 sin2 ( 2 </>))') 


(49 ) 

This form allows us to make several interesting ob¬ 
servations^ The crucial sign is that of 0 2 - In previous 
studiel^ES on iy 0 2 > 0 (so called phase A) was consid¬ 
ered. In our case however 0 2 < 0 and different ground 
state configurations should be considered. In phase B 
the minimization gives, n = ±m. Note two different so¬ 
lutions. So that the ’’vacuum manifold” is 


A = A 0 ne* x . 


(50) 


Here the range of x was enlarged, —n/2 < x < vr/2, 
to incorporate n = ±m. One can use the spherical rep¬ 
resentation of the unit vector (that will be termed ’’di¬ 
rector”): n = (sind n cos(p n ,sin0 n sin^„, cos0„). The 
ground state energy density therefore is achieved at 


a (T c - T ) 


0 1 + 02 


(51) 


Mathematically the vacuum manifold in phase B is iso¬ 
morphic to S 2 0 S 1 /Z 2 ■ This determines the thermody¬ 
namics of the superconductor very much in analogy with 
the scalar superconductor with 0 = + /? 2 . Magnetic 

properties are however markedly different. 


B. Small fluctuations analysis: two vastly different 
penetration depths and coherence lengths. 


Here the response of the superconductor in phase B 
to an external perturbation, like boundary or magnetic 
field, is considered. The basic excitation modes are un¬ 
covered by the linear stability analysis very similar to 
the so-called Anderson - Higgs mechanism in field theory 
applied to (scalar order parameter) superconductivity a 
long time ago 23 . Two basic scales, the coherence length 
(scale of variations of the order parameter) and mag¬ 
netic penetration depth (scale of variations of the mag¬ 
netic field), are obtained from the expansion of the GL 
energy to second order in fluctuations around supercon¬ 
ducting ground state at zero field. The fluctuations are 
parametrized by 6 n ,tp n ,Q m ,ip m ,<l>,Ai, A = A 0 (l+£). 


The order parameter, Eq.( 48), to the second order is 


A/A o «(O,O,l) + (0 n ,O,E + t$ (52) 

+ [ed n + OnVn - 0 2 / 2 + ie<p - </ 2 /2()53) 


A 4 


k 2 ( 0* n Q n + £*£ + (j>*cj) + i {e*4> - 4>*e)) 

+ (e* — kjAj + kjA* (e + i(j>) + A*Aj ; 


(55) 

(56) 


w = k\r n e n + k x k^i ( e + i<t>) + hk 3 (£ 

+k 2 ( e *£ + (j)*(f> + i — </>*£)) 

{ff‘ n A 1 + (e*-i</>*)A 3 + A* 1 e 
+ 3 V +A%{e + ict>) 


* 


n 


#*) 9n 
(57) 


+ a;a 3 . 


Diagonalization of this quadratic form reveals three 
Goldstone modes (that do not affect the characteristic 
lengths) and the ’’massive” fields, £ and A (the only one 
contributing for k = 0) with different longitudinal and 
transversal characteristic lengths: 


£2 _ U T p 2 _ U L C 2 . 

2 a(T c -T)’ 4L ut' 

\ 2 _ k 2 C 2 (0 s + 0 2 ) 2 _ UT 2 

T 8 tt e* 2 u T a(T c ~T)N' L u L T 


(58) 


Here £ T is the coherence length along the directions per¬ 
pendicular to the vector n, while the one parallel to n is 
£l- 

Similarly for magnetic penetration depths A t,l- Our 
calculation in the previous Section for the Dirac semi¬ 
metal, see Eq. (47), demonstrate that both are quite dif¬ 
ferent since ulJut = 1/32 << 1. This is obviously 
of great importance for large magnetic field properties 
of such superconductors. The phenomenological conse¬ 
quences for vortex state are briefly discussed in Section 
VI. 


C. Strong magnetic fields: is there an upper 
critical field H c 2? 


In strong homogeneous magnetic field H (assumed to 
be directed along the z axis) superconductivity typically 
(but not always, see an example of the p -wave supercon¬ 
ductor that develops flux phased) disappears at certain 
critical value H c2 . This bifurcation point is determined 
within the GL framework by the lowest eigenvalue of the 
linearized GL equations. This is a n exact requirement 
of stability of the normal phasd^SH. The linearized GL 
equation Eq.([32|) reads: 


(a - u T V 2 ) Sn 




A, = 0, 


(59) 


where coefficients are in Eq.(47), and u = ul — ut- We 
use the Landau gauge, A x = H c2 y; A y = A z = 0. As¬ 
suming translation symmetry along the field direction, 
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d~Ai = 0 , the operators of the eigenvalue problem de¬ 
pend on x and y only. 

Since we have three components of the order param¬ 
eter, there are three eigenvalues. It is easily seen from 
Eq.(591 that the z- component of the order parameter A z 
parallel to the external field direction is independent of 
the other two, A x , A y , leading to the ordinary Abrikosov 
value: 


- u t V 2 A z = -aA z -> Hi = 


(60) 


where £ T is defined in Eq.(58|. To avoid confusion with 
customary notations for layered materials (like high T c 
cuprates), the material that is modelled here is isotropic 
and ’’parallel”, ’’perpendicular” and refer to the relative 
orientation of the magnetic held to the vector order pa¬ 
rameter rather than to a layer. The orientation of the or¬ 
der parameter in isotropic material considered here, due 
to degeneracy of the ground state, is determined by the 
external magnetic held as we exemplify next. 

The two remaining eigenvalues involving only the order 
parameter components A x and A y perpendicular to the 
held are obtained from diagonalizing the ’’Hamiltonian”: 


higher r. In any case the Abrikosov lattice is expected 
to be markedly different from the conventional one and 
even from the vector order parameter model studied in 15 . 
At large helds in principle there is a direct interaction 
between the external magnetic held and spin, not taken 
into account in Hamiltonian Eq.Q studied so far. The 
next subsection addresses this question. 


D. What impact has the Pauli paramagnetism at 
strong fields? 

Since the prediction of the FFLO effect!^ in low T c 
superconductors it is well known that at very high mag¬ 
netic helds the direct spin - magnetic held coupling on 
the microscopic level might not be negligible. The sin¬ 
glet channel Cooper pair is effectively ’’broken” by the 
splitting since the spins of the two electrons are opposite 
(Pauli paramagnetic limit). It is not clear what impact it 
has on Dirac semi-metals. If the impact is large it could 
be incorporated as an additional paramagnetic term in 
the GL energy. In an isotropic Dirac superconductor one 
has only one possible term in the GL energy term linear in 
paramagnetic coupling and consistent with symmetries: 


n 



= —a 


A x 

A,, 


(61) 


(u T Vl + u L Vl %{D X ,T> V } \ 

{ %{D x ,Vy} u T Vl+u L VlJ- 


This nontrivial eigenvalue problem fortunately can be 
solved exactly, see Appendix C. The lowest eigenstate be¬ 
ing a superposition of just two lowest even Landau levels, 
|0) and \2) are: A x = a |0) + (3 |2), A y = 7 |0) + i/3 |2). 
The lowest of these eigenvalues is 


p* n^ 

e n c2 

he 


I ( u T +U L )_ 

— y/3 (Uj, + u 2 l ) - 2 U T U L 


= a (T c -T). 


(62) 

The corresponding critical held H^ (’’perpendicular” 
refers to the order parameter direction) that can be 
expressed via an effective ’’perpendicular” coherence 
length, 


rr±_ . 

11 c 2 — n .2 > 


2<r 
& = |(d + £)- 


(63) 

\/3Ci + (64) 


It is always larger than H'^ 2 , and therefore is physically 
realized. The upper held H^ becomes infinite at r c = 
u L /u T = (13 — 4\/l0) /3 ~ 0.117. This means that in 
such material superconductivity persists at any magnetic 
held like in some p- wave superconductors. It was found 
in Section IV that for the simplest Dirac semi-metal, r = 
1/32 < r c see Eq.(47). Thus there is no upper critical 
held in this case. Of course, different microscopic models 
that belong to the same universality class, might have 


F par = N^ p Ji( A*xA)-B, (65) 

where n p is the effective ” spin” of the Cooper pair some¬ 
times called ’’Zeeman coupling” ! 15 ! 24 !. 

The single particle Hamiltonian in magnetic held with 
the Pauli term becomes 

K = —ivFhV-a — /i + • B, ( 66 ) 

where the Bohr magneton, p B = eh/2mc, determines the 
strength of the coupling of the spin to magnetic held, with 
m being the free electron mass. The direct calculation, 
see Appendix A, shows that y p = 0. 


VII. DISCUSSION AND CONCLUSIONS 

Starting from the microscopic model of the isotropic 
Dirac semi-metal, the Ginzburg-Landau energy for this 
held is derived using the Gor’kov technique. It was found 
that the transversal coherence length is much smaller 
than the longitudinal, = \\/2^ T , despite the isotropy. 
Several new features appear when an external held is 
applied. The Ginzburg - Landau model is used to in¬ 
vestigate magnetic properties of such superconductors. 
Using the small deviation method the magnetic pene¬ 
tration depth was found also to be vastly different for 
longitudinal and transverse huctuations Xt/^l = 4-\/2. 

We have observed that properties of the triplet super¬ 
conductor phase of the Dirac senri-nretal has extremely 
unusual features that we would like to associate qualita¬ 
tively with the characteristics of the Cooper pair. The 












superconducting state generally is a Bose - Einstein con¬ 
densate of composite bosons - Cooper pairs, classically 
described by the Ginzburg - Landau energy as a func¬ 
tional of the order parameter. In the present case the 
Cooper boson is described by a vector field A,; (r). In this 
respect it is reminiscent to phonon and vector mesons in 
particle physical. Vector fields generally have both the 
orbital and internal degrees of freedom often called po¬ 
larization. The internal degree of freedom might be con¬ 
nected to the ’’valley” degree of freedom of constituents 
of the composite boson. We have provided evidence that 
the Cooper pair in DSM has finite orbital momentum, al¬ 
beit, as will be shown shortly, the spin magnetic moment 
is zero. Microscopically the unusual nature is related 
to the presence of the valley degeneracy in Dirac semi¬ 
metal. While in a single band superconductor the Pauli 
principle requires a triplet Cooper pair to have both odd 
angular momentum and spin, it is no longer the case in 
the Dirac semi-metal. 

A massive bosonic vector field in isotropic situa¬ 
tion (the case considered here) generally have distinct 
transversal and longitudinal polarizations (massless fields 
like photons in dielectric do not possess the longitudinal 
degree of freedom). The results for coherence lengths 

l and the penetration depths A t,l hr triplet super¬ 
conductor in DSM demonstrate pronounced disparity be¬ 
tween properties of transverse and longitudinal polariza¬ 
tions. In conventional ’’scalar” order parameter super¬ 
conductor the Abrikosov ratio®, n = A/£, distinguishes 
between type I (k < k c = 1 /V%) and type II, where, 
for example, vortices appear under a magnetic field. In 
the present case this separation is ambiguous. There 
are two quite different ratios: transversal and the lon¬ 
gitudinal, (in respect to director of the order parameter 
n = A/ |A|).see Eq.(47|. 


kt = - = 


e*UT 


he 




kl = 



(67) 

( 68 ) 


For typical DSM one estimates the Fermi velocity and 
chemical potential vp = c/200, p = 0.2eV, and with 
the expected critical temperatur^f® T c = 5A', one ob¬ 
tains the transversal Abrikosov ratio kt = 0.08, smaller 
than critical, while the longitudinal, up = 2.7 is larger. 
This means that longitudinal fluctuations the material 
behaves as type II, while response to the transverse ones 
is that of a type I superconductor. This has an impact 
on transport, optical and magnetic properties of these 
superconductors. 

The vortex physics of strongly type II triplet super¬ 
conductors of this type is very rich and some of it has al¬ 
ready been investigated in connection with heavy fermion 
and other superconductors suspected to possess p-wave 
pairing. In particular, their magnetic vortices appear as 


either vector vortices or so-called skyrmionei® - coreless 
topologically nontrivial textures. The magnetic proper¬ 
ties like the magnetization are very peculiar and even 
without a magnetic field the system forms a ’’sponta¬ 
neous flux state”. The material therefore can be called 
a ’’ferromagnetic superconductor”. The superconducting 
state develops weak ferromagnetism and a system of al¬ 
ternating magnetic domain^®. It was noted® that the 
phase is reminiscent to the phase B of superfluid He 3 !® 
(with an obvious distinction that the order parameter in 
the later case is neutral rather than charged and tenso- 
rial rather than vectorial). As was shown in section V, 
the Dirac semi-metal triplet superconductor phase is dif¬ 
ferent in several respects. It is more like the phase B of 
superfluid He 3 . 

Experimentally the major consequence of the present 
theoretical investigation, namely the polarization effect 
of the vector order parameter should be pronounced in 
the AC response of these materials. Recently the AC 
response of the disordered superconductor was utilized to 
probe Goldstone modeii®. We have demonstrated that 
they are abundant in the triplet DSM superconductor. 


VIII. APPENDIX A. 


A. Critical temperature calculation 


Starting from equation Eq. (26) the angle integrations 
result in (for p » Tp,T c ) 


1 _ T _ h 2 + _ 

9 ^ “n + {P 2 ~ h 2 f + Zu 2 n {p 2 + p 2 ) 

p 2 [ Td tanh (s/2T) _ p^_ 2 T d1e 

127r 2 J e =_t d £ 67r 2 ° S kTc 

where e = vfp — p. See the last (BCS) integral in 19 . 


(69) 


B. Cubic terms coefficients calculation 


To fix the two coefficients, /3 X and Eq. (32) we use 
only two components. The particular case j = k = l = 1 
(the coefficient of (/q 2 (/q) gives after angle integration 


nm /»00 

<™> 

p 2 (p 4 + 10p 2 (w 2 - 5 p 2 ) - 15 (p 2 + w 2 ) 2 ) 

—------2-wl) 

(p 4 + 2p 2 (w 2 - p 2 ) + (p 2 + w 2 ) 2 ) 

Performing finite integration (the upper bound on mo¬ 
mentum, p + Tp. can be replaced by infinity), one ob¬ 
tains 

Pi +^ 2 = Y5^4 S3> 
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where the sum is 


where 


1 7C (3) 


53 S(2« + l) 3 4 


(73) 


Similarly taking j = l = 2, k = 1 (the coefficient of 
V^Vh) gives after the angle integration 




/ p—0 


(74) 


P 2 (E 4 + 10p 2 (w 2 + 3p 2 ) + 15 (p 2 + w 2 ) 2 ) 
—- 2 —H75) 


(p 4 + 2p 2 (w 2 - /X 2 ) + (p 2 + W 2 ) 2 ) 
4p 2 


15?r 4 T 2 


S3- 


(76) 


resulting in Eq. (39). 


C. Effect of the Pauli interaction 


D" = 

*7 


( 82 ) 


(q 2 - {iuj + p) 2 ) 

{qjCii + qiOtj + Sij (ioj + n + a ■ q) + 2q z qjD}&3) 

In particular 

N (p) u L = - E Tr {T Z D (w, p) T*D" Z (w, p 

Ccip 




157r 4 T 2 uj??i 


S3 


7C (3) p 2 7£ (3) (vfK) N (p) 


607t 4 T 2 vf^i 


120 tt 2 T 2 


and 


TV 


(p) '«t = -^rr{r x D(a;,p)T*£)'' z (w,p)| = 32m l . 


cup 


The single particle Hamiltonian in magnetic held was 
written in Eq. (66). In order to fix the coefficient of the 


paramagnetic term linear in both the order parameter 
and Pauli coupling it is enough to expand the linearized 
Gorkov equations Eq. (35) to the first order in the spin 


density. Normal Greens functions have the following cor¬ 
rections: 


Dz ~ D — D (S • B) D ; 
D z « D + p s 5 (£* • B) D. 


(77) 


The Pauli term in Gor’kov equation (aft er m ultiplying 
by T\ and taking the trace as usual), Eq.(32), therefore 
is obtained from expansion of Eq. (35), 


^Tr|T z DA*£>| 


(78) 


CO p 


= -iv z e ijk A*B k N (p) = » B B? jk A*B k , (79) 

Bf jk = E Tr { T ' b - T i D ^) d } - 


The bubble sum is directly evaluated and vanishes Bf ;jk = 
0 . 


X. APPENDIX C. EXACT SOLUTION FOR 
UPPER CRITICAL MAGNETIC FIELD 


In this Appendix the matrix T~L defined in Eq.(61) de¬ 


termining the perpendicular upper critical held is diago¬ 
nalized variationally. 


A. Creation and annihilation operators 

Using Landau creation and annihilation operators in 
units of magnetic length = l~ 2 for the state with 
k x = 0 (independent of a:), so that covariant derivatives 
are 


D x = d x +iy = iy = -j=(a + a + ) ; (84) 


'y — Vy 


1 

3 


In terms of these operators the matrix operator 7~L takes 
a form: 


IX. APPENDIX B. CALCULATION OF 
GRADIENT TERMS 

Rotational invariance allows to represent the sum in 


Eq (46) terms of coefficients ut and ul- 

(p) (ut ( P 2 5mj — PmPj ) + U L PmPj) (80) 


= PkPl 


E Tr ^T f m D (lv. q) Tj*D kl (uj, q 

cuq 


(81) 


n=u T + - + V; (85) 

Vn = 2 ura + a + ^ (a 2 + a +2 + 2a + a) (86) 


tu 


V 12 = V 21 = f (a+ 2 - a 2 ) 


(87) 


V 22 = 2 w T a + a — — (a 2 + a +2 — 2a + a) (88) 

The exact lowest eigenvalue is a combination of two 
lowest Landau levels. Indeed applying the operator V on 
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a general vector on the subspace gives 


eigenvalue equations 


V 


a\0}+(3\2}\ 

7 |0) + <512) ) 


/ P) | 0 ) \ 

-H^ + S +2/3 ) +4 ^) l 2 > 

+ f (~P _ i$) 14) 

+^(i/? + <5)|0> 

+ ( u (-^ + ^- 2<5 )-M |2) 

V +f(-t/3 + 5)|4) / 


(89) 


For S = i/3, higher Landau levels decouple and one gets 


—v 


V2 

0 


— y/2u 

—Aut — 2 u — v 

uiV 2 


0 

_ in 
\/2 
— V 


= 0 , 


resulting in three eigenvalues of Lf 

M 1 ) = «t + u/2, (90) 

= 3 ut + —u ± \J 4+ 4 utu + 3 u 2 . (91) 
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